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Positive  polarity  applied  voltage  and  gas  pressure  dependent  scaling  of  cathode  directed  streamer 
propagation  properties  in  helium  gas  flow  guided  capillary  dielectric  barrier  discharge  have  been 
quantifled  from  streamer  velocity,  streamer  current,  and  streamer  optical  diameter  measurements. 

All  measurements  of  the  non- stochastic  streamer  properties  have  been  performed  in  a  variable 
gas  pressure  glass  cell  with  N2  co-flow  and  under  self-consistent  Poisson  electric  field  dominated 
conditions  to  permit  data  comparison  with  2-D  streamer  dynamics  models  in  air/nitrogen.  The 
streamer  optical  diameter  was  found  to  be  nearly  independent  of  both  gas  pressures,  from  ITOTorr 
up  to  760Torr,  and  also  for  applied  voltages  from  6  to  llkV  at  520Torr.  The  streamer  velocity 
was  found  to  increase  quadratically  with  increased  applied  voltage.  These  observed  differences  in 
the  2-D  scaling  properties  of  ionization  wave  sustained  cathode  directed  streamer  propagation  in 
helium  flow  channel  with  N2  annular  co-flow  compared  to  the  streamer  propagation  in  air  or  nitro¬ 
gen  have  been  shown  to  be  caused  by  the  remnant  ionization  distribution  due  to  large  differences  in 
the  dissociative  recombination  rates  of  He2^  versus  N4^  ions,  for  this  5  kHz  repetition  rate  applied 
voltage  pulse  generated  streamers.  ©2014  AIP  Publishing  LLC. 

[http://dx.doi.Org/10.1063/l.4892968] 


1.  INTRODUCTION 

Highly  non-equilibrium  atmospheric  pressure  plasma 
jets  have  recently  been  investigated  for  producing  active 
plasmas  under  ambient  conditions.  One  specific  type  of  such 
plasma  jet  that  allows  the  formation  of  a  non- stochastic 
active  discharge  extending  beyond  the  anodic  electrode  uses 
a  capillary  dielectric  tube  with  helium  gas  flow  under  ambi¬ 
ent  conditions.  Several  experimental^”^  and  modeling^”^"^ 
results  of  such  helium  plasma  jets  have  been  published 
which  show  that  this  type  of  discharge  behaves  like  a  cath¬ 
ode  directed  streamer  propagation.  This  plasma  jet  uses  a 
dielectric  barrier  configuration  with  annular  anode  config¬ 
ured  as  a  co-axial  electrode  with  the  helium  gas  flow,  in  con¬ 
trast  to  a  point-to-plane  electrode  configuration  normally 
used  for  streamer  propagation  in  air  or  in  nitrogen.  The  gas 
dielectric  breakdown  and  self- sustained  streamer  propaga¬ 
tion  are  possible  to  achieve  in  the  helium  gas  channel  at  a 
much  lower  reduced  electric  field  compared  to  the  co¬ 
flowing  and/or  surrounding  air.  These  differences  in  electri¬ 
cal  characteristics^^  of  streamer  propagation  in  helium  versus 
air^^’^^  are  responsible  for  the  observed  differences  in  the 
propagation  of  ionization  wave  sustained  streamer  propaga¬ 
tion  in  helium  flow  channel  in  ambient  air  compared  to  the 
streamer  propagation  in  air  or  nitrogen.  Although  many 
experimental  results  and  model  calculations  have  been  pub¬ 
lished  for  streamer  propagation  in  the  helium  flow  channel  in 
ambient  air,  a  systematic  study  of  applied  voltage  and/or 
pressure  dependent  two-dimensional  or  three  dimensional 
scaling  properties  of  the  streamer  launched  by  this  capillary 
dielectric  barrier  discharge  (CDBD)  with  annular  anodic 
electrode,  compared  to  point-to-plane  streamer  propagation 


in  air,  have  not  yet  been  published.  In  this  work,  we  report 
measurements  of  applied  voltage  and  gas  pressure  dependent 
scaling  properties  of  streamer  current,  streamer  speed,  and 
streamer  radiative  channel  diameter with  co-flowing 
nitrogen. 

In  order  to  control  the  purity  of  the  co-flowing  gas  and 
avoid  the  effects  of  ambient  room  air  current,  the  capillary 
dielectric  tube  with  annular  anode  was  placed  inside  a  3  cm 
diameter  glass  enclosure  where  the  gas  pressure  was  con¬ 
trolled  by  a  downstream  pump  and  flow  controller.  One  of 
the  advantages  of  controlling  the  co-flow  gas  is  that  the 
plasma  jet  propagation  can  be  maintained  over  a  relatively 
long  length  as  the  total  gas  pressure  is  varied  from  760Torr 
down  to  170Torr.  This  glass  enclosed  cell  experimental  con¬ 
dition  is  very  well  suited  for  the  systematic  study  of  the 
streamer  scaling  properties  in  lower  dielectric  breakdown 
strength  gases  surrounded  by  much  higher  dielectric  break¬ 
down  strength  gases.  The  ring  anode  electrode  used  in  this 
CDBD  configuration  produces  a  Laplacian  electric  field  in 
addition  to  the  self-consistent  Poisson  electric  field  of  the 
streamer.  In  order  to  study  the  voltage  and  gas  pressure  de¬ 
pendent  scaling  properties  of  streamer  propagation  under  the 
influence  of  low  external  applied  electric  field  and  streamer 
dynamics  properties  dominated  by  the  self-consistent 
Poisson  electric  field,  it  is  required  to  measure  streamer  scal¬ 
ing  properties  at  locations  far  away  from  the  ring  anode. 
This  glass  cell  enclosed  plasma  jet  configuration  permits  us 
to  measure  streamer  velocity,  streamer  radiative  channel  di¬ 
ameter,  and  streamer  current  scaling  with  the  variation  of 
both  applied  voltage  and  gas  pressure  under  the  conditions 
of  low  applied  external  electric  field.  We  have  used  spatio- 
temporally  resolved  imaging  of  spectrally  resolved  plasma 
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emission  from  the  helium  3^D— 2^P  transition  to  measure 
streamer  velocity  and  streamer  radiative  channel  diameter. 
We  have  also  measured  spatially  (1  mm  at  the  focal  volume) 
and  temporally  (±  2  ns)  resolved  optical  emission  from  He 
3^S-  2^P  at  706.5  nm,  391  nm  N2+  and 

337.1  nm  emission  from  N2  C^Xlu-B^rig  using  optical  fiber 
coupled  photomultipliers.  Of  these,  the  587.6  and  706.5  nm 
emission  lines  from  atomic  helium  provided  accurate  signa¬ 
tures  for  direct  electron  impact  ionization  of  helium  at  the 
streamer  head.  If  the  direct  electron  impact  ionization  of  he¬ 
lium  and  nitrogen  gases  are  occurring  at  the  streamer  head,^^ 
then  dual  optical  hber  channel  coupled  photomultipliers 
should  detect  those  signatures  from  our  simultaneous  optical 
emission  intensity  measurements. 

Experimental  details  and  experimental  results  of  two- 
dimensional  scaling  behavior  of  the  streamers  in  lower 
dielectric  breakdown  strength  helium  gas  surrounded  by 
much  higher  dielectric  breakdown  strength  nitrogen  gas  are 
discussed  below.  A  zero  dimensional  plasma  chemical 
kinetics  model  has  been  used  to  show  the  influence  of  rem¬ 
nant  ionization  on  the  streamer  propagation  properties  of  the 
non- stochastic  helium  plasma  jet  for  5  kHz  repetition  rate 
pulsed  excitation. 

II.  EXPERIMENTAL  DETAILS 

Figure  1  shows  the  helium  plasma  jet,  powered  by  a  sin¬ 
gle  copper  ring  anode,  mounted  coaxially  inside  a  30  mm 
inner  diameter  cylindrical  glass  cell  having  a  co-flow  of 
nitrogen.  Downstream  from  the  plasma  jet  assembly,  the  cell 
gas  pressure  was  automatically  monitored  and  maintained  at 
pressures  ranging  from  760Torr  down  to  170Torr  while  a 


FIG.  1.  Visible  emission  photograph  of  the  2  SLM  helium  plasma  jet  operat¬ 
ing  within  the  glass  cell  with  a  0.5  SLM  N2  annular  co-flow  at  520Torr 
excited  by  9  kV  unipolar  total  applied  voltage  pulses. 


mechanical  vacuum  pump  was  used  to  evacuate  the  gases. 
The  glass  capillary  was  6  cm  long  with  an  inner  diameter  of 
2  mm  and  an  outer  diameter  of  3  mm.  As  shown  in  Figure  1 , 
a  2  cm  long  copper  ring  electrode  was  placed  2  mm  from  the 
capillary  tip.  For  all  conditions  reported  in  this  paper,  helium 
gas  flowed  through  the  capillary  at  2  standard  liters  per  mi¬ 
nute  (SFM)  resulting  in  an  exit  flow  speed  of  ^11  m/s.  The 
co-flowing  nitrogen  background  gas  was  introduced  into  the 
cell  25  cm  behind  the  capillary  tip  at  a  volume  flow  rate  of 
0.5  SFM.  Based  on  the  inner  diameter,  the  Reynolds  number 
for  the  helium  flow  at  the  capillary  tube  exit  was  about  190 
at  520  Torr. 

Unipolar,  high  voltage  pulses  (6-12  kV)  having  a  20  ns 
rise  time  were  applied  to  the  copper  ring  anode  at  a  5  kHz 
pulse  repetition  rate  via  an  electrical  feed  through  in  the 
glass  cell  wall.  For  the  time-averaged  streamer  plasma  emis¬ 
sion  image  at  520  Torr,  shown  in  Figure  1,  the  peak  total 
applied  voltage  at  the  anode  was  9  kV  measured  at  the  high 
voltage  switch,  prior  to  a  300  Q  ballast  resistor,  using  a 
1000:1  high  voltage  probe  (North  Star  PVM-11).  The  total 
circuit  current  was  measured  at  the  powered  anode  wire 
using  an  annular  wide-bandwidth  current  sensor  (Pearson 
model  4100,  1  V/A).  Also  shown  in  Figure  1,  the  streamer 
current  was  measured  by  symmetrically  threading  the  glass 
cell  through  the  5  cm  inner  diameter  hole  of  an  annular  wide- 
bandwidth  current  sensor  (Pearson  model  6585,  1  V/A). 
Unless  otherwise  noted,  the  front  face  of  the  streamer  current 
sensor  was  located  23  mm  from  the  capillary  tip  for  all  ex¬ 
perimental  conditions.  The  time- averaged  broad  band  visible 
emission  from  the  plasma  jet,  as  shown  in  Figure  1,  implies 
that  the  streamer  discharge  propagated  through  the  25  mm 
wide  current  probe  and  then  at  least  another  2  cm  further 
downstream  under  this  gas  flow  condition.  All  time- 
dependent  voltage  and  current  measurements  were  acquired 
with  a  4-channel,  500  MHz  bandwidth  oscilloscope  (FeCroy 
6050),  1000  samples  averaged,  and  then  transferred  to  a 
desktop  computer  for  analysis. 

A  gated  intensified  charge-coupled  device  (ICCD)  cam¬ 
era  (Princeton  Instruments  PI-MAX  I,  512x512  CCD  pixel 
arrays)  was  used  to  collect  spatiotemporally  resolved  atomic 
helium  3^D-2^P  emission  at  587.6  nm  using  a  10  nm  band 
pass  filter  (Thorlabs  FB590-10-1)  mounted  at  the  front  of  the 
camera  lens  system.  The  He  3^D— 2^P  transition  was  chosen 
because  it  has  a  direct  electron  impact  excitation  threshold 
energy  of  23.1  eV,  compared  to  24.6  eV  for  ionization,  and  a 
relatively  large  spontaneous  emission  coefficient  (Ai2~ 
7  X  10^s“\  Ref.  27),  making  it  a  sensitive  marker  for  direct 
electron  impact  excitation  at/or  near  the  streamer  head.  The 
imaging  optical  axis  was  perpendicular  to  the  capillary  axis. 
The  camera  collection  optics  consisted  of  a  BK-7  telecentric 
and  zoom  lens  combination  with  several  extender  rings  and  a 
working  distance  of  170  mm,  resulting  in  a  1:1  magnification 
of  the  optical  collection  system.  The  12  mm  field  of  view 
along  the  streamer  propagation  axis  was  between  6  mm  and 
18  mm  from  the  capillary  tip,  and  the  image  resolution  was 
^150 /im,  corresponding  to  ^6  pixels.  Unless  otherwise 
noted,  all  images  were  acquired  using  a  5  ns  ICCD  gate  on- 
time,  40000  accumulated  ICCD  exposures,  and  a  fixed 
intensifier  gain  setting. 
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III.  RESULTS 
A.  520  Torr  data 

We  have  chosen  to  use  520  Torr  closed  cell  pressure  for 
initial  sets  of  measurements  since  it  allowed  us  to  measure 
streamer  scaling  properties  over  a  wide  voltage  range  with¬ 
out  surface  flashover  in  our  conhguration.  Compared  to  the 
open  plasma  jet  configuration  at  760  Torr,  it  also  allowed  us 
to  control  the  co-flowing  gas  compositions. 

Figure  2  shows  a  typical  example  of  the  current  and 
voltage  traces  acquired  for  the  plasma  jet  operating  at  a  peak 
total  applied  voltage  of  9kV  using  N2  as  the  co-flowing 
background  gas.  Note  that  the  timescale  in  Figure  2  is  rela¬ 
tive  to  the  peak  total  circuit  displacement  current.  The  total 
circuit  current  and  applied  voltage  were  synchronized  to  an 
accuracy  of  ±2  ns  by  temporally  shifting  and  rescaling  the 
time  integrated  displacement  current  and  matching  the  am¬ 
plitude  to  that  of  the  applied  voltage  temporal  signal.  The 
streamer  current  was  obtained  by  first  acquiring  a  back¬ 
ground  trace  from  the  streamer  current  probe  with  the  helium 
flow  turned  off,  but  with  the  applied  voltage  on  (discharge 
off  condition).  This  background  trace  was  stored  in  the  oscil¬ 
loscope  memory  and  then  subtracted  from  the  real  time 
streamer  probe  signal.  With  the  resumption  of  helium  flow, 
the  discharge  was  turned  back  on  and  steady  state  conditions 
were  established.  Ten  real  time  streamer  current  samples 
were  acquired  and  then  averaged  during  post  processing. 
Non-random  noise  fluctuations  “ringing”  due  to  circuit 
effects  were  mitigated  using  adjacent  average  time-step  fil¬ 
tering  (OriginPro  8.1  software  package^^).  For  each  dis¬ 
charge  condition,  the  signal  filtering  was  carefully  chosen  to 
preserve  both  the  rise  time,  peak,  and  decay  time  features. 
The  signal  propagation  delay  with  respect  to  the  total  circuit 
current  was  accounted  for  by  temporally  shifting  the  onset 
time  of  the  streamer  probe  background  trace  so  that  it 
matched  the  onset  time  of  the  total  circuit  displacement  cur¬ 
rent  with  about  ±3  ns  accuracy.  This  delay  time  correction 
was  applied  to  the  mean  streamer  current  trace.  For  the  9kV 
peak  applied  voltage,  in  Figure  2,  the  measured  peak 
streamer  current  was  20  mA,  and  an  estimated  upper  bound 
of  60  jiiJ  energy  was  deposited  into  the  streamer. 


10n 


8- 


> 

0) 

O) 

CD 

?  6 
> 

■D 

4 

Q. 

Q. 

i  2^ 

o 

0-E 


-Total  Applied  Voltage 
-Total  Circuit  Current ^  200 
-  Streamer  Current 


i25 


20 


15  I 

10  i 

3 

o 

-  5 


0.0 


0.4  0.8 

Time  (^s) 


1.2 


FIG.  2.  Measured  total  applied  voltage  (dashed  line),  total  circuit  current, 
and  streamer  current  (dotted  line)  for  the  helium  plasma  jet  excited  by  9  kV 
peak  unipolar  high  voltage  pulses  at  5  kHz  pulsed  repetition  rate. 


Streamers  generated  by  CDBD  helium  plasma  jets 
operating  in  air  are  known  to  exhibit  an  annular  morphol¬ 
ogy  having  an  outer  diameter  that  is  initially  approximately 
equal  to  the  capillary  tube  inner  diameter. The  ring¬ 
like  shape  is  formed  due  to  the  initial  Laplacian  electric 
field  that  induces  breakdown  near  the  inner  wall  of  the 
dielectric  tube.^’^^’^"^’^^  As  the  streamer  propagates  away 
from  the  capillary  tip,  confined  within  the  helium  gas  jet, 
the  annular  radius  tends  to  decrease  towards  the  centerline 
axis  and  may  eventually  converge  on  axis.^’^^’^^’^^  From 
the  two  dimensional  ICCD  emission  images,  radial  inten¬ 
sity  slice  plots  were  obtained  using  Image  J  imaging  analy¬ 
sis  software^^  for  applied  voltages  from  6-10  kV.  Both  the 
two  dimensional  emission  intensity  images  and  radial  slice 
plots  showed  that  the  streamer  head  exhibited  a  symmetric 
annular  structure  that  converged  on-axis  approximately 
8  mm  from  the  capillary  tip.  Axial  slice  plots  obtained  from 
late  time  emission  intensity  images  (not  shown)  for  which 
the  streamer  had  moved  well  past  the  field  of  view  showed 
on-axis  He  587.6  nm  emission  attributed  to  the  Laplacian 
electric  field  extending  beyond  8  mm  up  to  approximately 
12  mm  (recall  that  the  total  applied  voltage  fall  time  was 
^8  fis).  Taken  together,  the  information  obtained  from 
these  ICCD  images  indicates  that  the  streamer  head  radial 
emission  profile  exhibited  a  morphological  transition  from 
ring-like  near  the  capillary  tip  to  Gaussian-like  at  16  mm 
from  the  capillary  tip.  At  this  16  mm  location,  the  influence 
of  the  Laplacian  electric  field  on  the  streamer  head  was 
negligible  (discussed  below  in  Sec.  Ill  B)  and  behaves 
more  like  a  cathode  directed  positive  streamer  with  low 
applied  electric  field  (in  Sec.  Ill  B,  we  illustrate  this  transi¬ 
tion  in  Figure  8  which  shows  ICCD  radial  emission  inten¬ 
sity  slice  plots  at  6  mm  and  16  mm  axial  positions  obtained 
at  365  Torr  and  6kV  applied  voltage).  For  this  reason,  we 
have  chosen  the  16  mm  axial  location  for  comparing 
streamer  scaling  properties  with  model  calculations  in  he¬ 
lium  and  and  also  measurements in  air 

and  nitrogen,  where  streamer  scaling  properties  are  deter¬ 
mined  by  space-charge  induced  Poisson  electric  field  at  the 
streamer  head. 

Figure  3  shows  the  scaling  of  peak  streamer  current  with 
peak  total  applied  voltage  taken  at  two  different  locations  of 
the  current  probe  at  520  Torr.  The  square  symbols  represent 
the  data  with  the  current  sensor  front  face  located  23  mm 
from  the  capillary  tip.  The  circles  represent  the  peak  current 
scaling  for  conditions  with  the  front  face  of  the  Pearson 
probe  located  33  mm  from  the  capillary  tip.  Also  shown  in 
Figure  3  are  best-fit  polynomial  lines  through  the  data  points 
which  indicate  a  quadratic  scaling  of  peak  current  increase 
with  increasing  voltage.  We  define  the  meaning  of  the  “best 
fit”  as  that  polynomial  expression  which,  statistically,  best 
describes  the  data  set  using  the  lowest  degree  possible  with  a 
95%  confidence  level  constraint.  The  assessment  of  the  mini¬ 
mal  degree  required  was  performed  by  sequentially  testing 
four  sets  of  polynomials  ranging  from  fourth  order  down  to 
first  order  using  Analysis  of  Variance  (ANOVA)  error- 
weighted  polynomial  regression.^^’^^  The  quadratic  depend¬ 
ence  of  streamer  current  with  applied  voltage  shown  in 
Figure  3  is  in  agreement  with  1-D  analysis  and  2-D  modeling 
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FIG.  3.  Measured  peak  streamer  current  as  a  function  of  total  applied  volt¬ 
age  at  520Torr  with  the  probe  front  face  located  23  mm  (squares)  and 
33  mm  (circles)  from  the  capillary  tip.  The  solid  lines  are  quadratic  best  fits 
(see  text)  to  the  data  points. 

for  positive  streamer  propagation  in  homogeneous  air/^  As 
the  current  probe  was  moved  from  23  mm  to  33  mm  from  the 
capillary  tip,  the  streamer  current  remains  nearly  unchanged. 
This  suggests  that  diffusional  mixing  of  the  annular  co-flow 
gas  into  the  helium  core  flow  can  be  neglected  at  positions 
approximately  up  to  33  mm  from  the  capillary  tip. 

Figure  4  shows  the  scaling  of  streamer  velocity  increase 
with  total  peak  applied  voltage  when  the  maximum  streamer 
head  emission  intensity  was  located  between  15  and  16  mm 
from  the  capillary  tip.  These  measurements  were  performed 
simultaneously  with  the  peak  streamer  current  measurements 
shown  in  Figure  3.  The  streamer  speeds  for  each  discharge 
condition  were  derived  from  sequentially  time  delayed  gated 
ICCD  camera  images  of  the  He  587.6  nm  emission  at  the 
streamer  head  within  the  imaging  held  of  view  between 
6  mm  and  18  mm  from  the  capillary  tip.  Using  a  5  ns  gate 
width,  the  intensifler  gate  delay  was  adjusted  so  that  at  least 
three  (depending  upon  the  streamer  speed)  background  sub¬ 
tracted  images  having  a  temporal  separation  no  more  than 
5  ns  apart  were  acquired.  For  this  2-D  axisymmetric  geome¬ 
try,  the  streamer  head  propagation  direction  is  collinear  with 
the  electric  held,  so  the  measured  streamer  speed  should  be 


same  as  the  streamer  head  velocity.  From  these  images,  10- 
pixel  radial  binned  average  intensity  slice  plots  cutting 
across  the  entire  image  along  the  propagation  axis  and 
through  the  streamer  head  axial  centerline  were  obtained 
using  ImageJ.  From  these  emission  intensity  slice  plots,  the 
average  streamer  velocities  were  estimated  by  dividing  the 
change  between  the  90%  peak  height  axial  positions  of  tem¬ 
porally  successive  slice-plot  pairs  by  the  change  in  their  re¬ 
spective  gate  delay  times.  It  should  be  noted  that  the 
estimated  streamer  velocities  at  each  applied  voltage  were 
nearly  constant  (within  the  10%  uncertainty  of  the  measure¬ 
ments)  between  14  mm  and  17  mm  from  the  capillary  tip.  As 
shown  in  Figure  4,  by  the  minimum  polynomial  fit  of  the 
data,  the  streamer  velocity  increase  scaled  quadratically  with 
increased  applied  voltage.  Both  measurements^^  and  model¬ 
ings^  in  air  have  predicted  a  linear  relationship  between  posi¬ 
tive  streamer  velocity  and  applied  voltage  for  streamer 
propagation  in  homogeneous  air.  It  is  not  known  if  our  mea¬ 
surement  results,  which  indicate  quadratic  streamer  velocity 
scaling  with  increased  applied  voltage  is  due  to  a  departure 
from  the  velocity-diameter  scaling  (see  next  paragraph) 
observed  for  streamer  propagation  in  air^s)^2i,26 
“streamer  coupling  for  glow  like  discharge’’^^  in  the  helium 
gas  channel  with  high  degree  of  remnant  ionization  in  the 
5  kHz  pulse  repetition  rate  condition  (see  Sec.  IV). 

Figure  5  shows  a  plot  of  the  streamer  head  radiative 
channel  “optical”  diameter  (FWHM)  located  16  mm  from 
the  capillary  tip  as  a  function  of  total  applied  voltage 
obtained  from  ICCD  image  emission  intensity  radial  slice 
plots  (radial  profiles  obtained  at  6kV  and  8kV  at  520Torr 
are  shown  in  Figures  9(a)  and  9(b)).  The  set  of  optical  diam¬ 
eter  measurement  values  are  also  listed  in  Table  I.  We  note 
that  these  measurements  were  performed  simultaneously 
with  the  streamer  current  measurements  shown  in  Figure  3 
with  the  probe  front  face  located  23  mm  from  the  capillary 
tip  and  were  also  obtained  from  the  same  images  used  to 
estimate  the  streamer  speeds  shown  in  Figure  4.  Both  Figure 
5  and  Table  I  show  that  the  streamer  head  optical  diameters 
are  essentially  independent  of  the  applied  voltage.  Note  that 
the  optical  diameters  change  14%  from  0.63  mm  at  6kV  to 
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FIG.  4.  Streamer  velocity  as  a  function  of  total  applied  voltage  at  520Torr. 
The  streamer  head  was  located  16  mm  from  the  capillary  tip.  The  solid  line 
is  a  quadratic  best  fit  to  the  data  points  (see  text). 


FIG.  5.  Streamer  optical  diameters  (FWHM)  as  a  function  of  total  applied 
voltage  at  520  Torr.  The  streamer  head  was  located  16  mm  from  the  capillary 
tip.  The  solid  line  is  a  linear  best  fit  to  the  data  points  (see  text). 
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TABLE  I.  Optical  diameters  at  the  streamer  head  measured  16  mm  from  the 
capillary  tip  using  ICCD  imaging  of  the  He  587.6 nm  emission  intensity. 
Also  shown  are  electron  number  density  estimates  as  a  function  of  applied 
voltage  at  520Torr.  The  uncertainty  of  the  diameter  measurements  are 
±13%  while  that  of  the  number  density  estimates  are  ±20%. 


Applied  voltage  (kV) 

D  (mm) 

He  (x  10^^  cm 

6 

0.63 

5.5 

7 

0.68 

6.3 

8 

0.73 

6.4 

9 

0.73 

6.7 

10 

0.73 

7.5 

0.73  mm  at  lOkV  while  the  corresponding  streamer  veloc¬ 
ities  shown  in  Figure  4  exhibit  a  70%  increase  from 
0.20  mm/ns  to  0.61  mm/ns  between  6  and  10  kV.  The  lower 
values  of  optical  diameter  at  6  kV  and  7  kV  applied  voltages 
(see  Table  I)  are  attributed  to  the  lower  ICCD  imaging  S/N 
ratio.  As  mentioned  above,  the  radial  ICCD  emission  inten¬ 
sity  profiles  were  Gaussian-like  at  the  16  mm  position  (for 
example,  see  Sec.  Ill  B  Figures  8,  9(a),  and  9(b)).  The  diam¬ 
eters  near  the  baseline  using  a  2  ns  ICCD  gate  were  approxi¬ 
mately  0.2  mm  smaller  than  those  acquired  using  a  5  ns  gate 
time,  a  difference  that  is  near  the  imaging  system  resolution 
and  can  be  attributed  to  a  lower  S/N  ratio  rather  than  the  effect 
of  “blurring”  due  to  axial  translation  during  the  exposure 
time.  The  FWHM  optical  diameters  were  only  sensitive  to  the 
2  ns  ICCD  gate  time  at  6  kV  and  7  kV  and  were  again  within 
the  imaging  system  resolution.  We  chose  the  FWHM  optical 
diameter  to  be  the  metric  for  the  streamer  head  diameter 
because  of  the  better  S/N  ratio  and  as  a  standard  for  compari¬ 
son  to  other  streamer  data  in  the  literature. 

Figure  6  shows  the  velocity-diameter  plot  for  the  same 
data  taken  from  Figures  4  and  5.  Because  the  optical  diame¬ 
ters  are  nearly  constant  as  a  function  of  applied  voltage  (see 
Table  I  and  Figure  5),  the  streamer  velocity  is  independent 
of  optical  diameter  for  the  applied  voltage  range  of  6-10  kV. 
This  velocity-diameter  scaling  in  the  helium  plasma  jet 
exhibits  a  strong  departure  from  measurements^"^’^^  and  mod¬ 
eling  results of  positive  streamers  in  air.  The  invariance  of 
the  optical  diameter  with  applied  voltage  and  streamer  ve¬ 
locity  is  a  consequence  of  both  the  lower  dielectric 
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FIG.  6.  Streamer  velocity  as  a  function  of  optical  diameter  at  520  Torr.  The 
streamer  head  was  located  16  mm  from  the  capillary  tip. 


breakdown  strength  of  the  helium  gas  channel  compared  to 
the  surrounding  nitrogen  co-flow,  which  has  a  much  higher 
breakdown  electric  field  requirement but,  more  impor¬ 
tantly,  it  may  be  due  to  remnant  ionization  in  pure  helium 
core  fiow  versus  helium-nitrogen  boundary  layer  for  this 
5  kHz  repetition  rate  pulsed  discharge. 

We  were  not  able  to  quantify  accurately  the  scaling  of 
the  streamer  head  emission  axial  FWHM  with  applied  volt¬ 
age  from  the  ICCD  images  due  to  axial  “blurring”  of  the 
propagating  streamer  head  within  the  5  ns  ICCD  camera  gate 
time.  However,  we  have  used  images  taken  with  5  ns  and 
2  ns  gate  times  and  extrapolated  the  axial  FWHM  to  near 
0  ns  gate  time  for  the  condition  when  the  current  probe  front 
face  was  located  23  mm  from  the  capillary  tip.  Assuming  a 
linear  relationship  between  the  axial  FWHM  and  gate  time, 
an  extrapolation  to  0  ns  gate  time  resulted  in  an  axial  FWHM 
between  0.6  mm  at  8kV  to  1.4  mm  at  10  kV.  Again,  these 
characteristics  of  the  radial  versus  axial  electron  impact  exci¬ 
tation  profiles  in  helium  seem  to  be  different  compared  to 
•  •  1 8—20 
streamers  m  air. 

Electron  number  density  estimates  behind  the  streamer 
head  for  each  applied  voltage  pulse  have  been  made  using 
the  one  dimensional  current  continuity  approximation. 


4/ 


rie  = 


evnD^  ’ 


(1) 


where  /,  v,  and  D  represent  the  measured  values  of  the 
streamer  currents,  streamer  velocities,  and  optical  diameters 
given  in  Figures  3,  4,  and  5,  respectively.  With  this  approxi¬ 
mation,  we  have  assumed  that  (i)  the  radial  component  of  the 
secondary  electron  drift  flux  is  not  significant,  (ii)  the 
streamer  velocity  is  equal  to  the  axial  electron  drift  velocity, 
and  (iii)  the  measured  current  is  uniformly  distributed  over 
the  area  of  the  optical  diameter.  Under  these  assumptions, 
estimates  of  the  electron  number  densities  vs.  applied  volt¬ 
age  using  Eq.  (1)  are  given  in  Table  I.  These  electron  number 
density  estimates  range  from  approximately  5  x  10“  to 
7  X  10“cm“^  for  increasing  applied  voltage  with  an  uncer- 
tainty  of  20%.  The  maximum  electric  field  within  the 
streamer  head,  E*,  can  be  related  to  the  channel  electron 
density  behind  the  streamer  head  through  the  expression  n^ 
^  (ao/e)a*E*,^^’^^  where  a*  is  the  ionization  coefficient  for 
helium  evaluated  for  Using  this  expression  with  the  val¬ 
ues  of  number  density  given  in  Table  I  along  with  the  func¬ 
tional  dependence  of  a  on  the  reduced  electric  field 
for  pure  helium,^^  we  estimate  that  E*  is  approximately 
10kV/cm-12kV/cm  at  520  Torr  (60-70  Td)  for  the  applied 
voltages  range  of  6-10  kV.  This  peak  electric  field  estimate 
is  in  agreement  with  modeling  results^^’^^  and  experiment^^ 
for  helium  gas  guided  streamers  at  atmospheric  pressure. 
Note  that  modeled  peak  electric  field  for  streamers  in  air  at 
atmospheric  pressure^ are  approximately  an  order  of 
magnitude  larger  with  corresponding  electron  densities  on 
the  order  of  10'"-10'^cm-l^'’^^-^^’^" 

We  have  also  measured  spatially  (1mm  at  the  focal 
volume)  and  temporally  (±2  ns)  resolved  optical  emission  from 
He  3^D— 2^P  transition  at  587.6  nm.  He  3^S-  2^P  at  706.5  nm, 
N2+  391  nm  (BX+- 


^  X  Eg^),  and  337.1  nm  emission  from 
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N2  C^rTu-B^rig  using  dual  fiber  coupled  photomultipliers. 
The  temporal  variation  of  He  706.5  nm,  N2^  391  nm,  and  N2 
337.1  nm  emissions  has  been  measured  with  respect  to  the 
He  587.6  nm  emission  line.  The  normalized  intensity  plots 
for  all  those  transitions  are  shown  in  Figure  7.  In  Figure  7, 
all  traces  have  been  corrected  for  cable  signal  propagation 
delays  relative  to  the  He  587.6  nm  emission  trace  but  not  for 
the  photo  multiplier  tube  (PMT)  transit  time  delays  (47  ns). 
A  correction  for  PMT  transit  time  delays  would  result  in  a 
47  ns  shift  of  all  waveforms  earlier  in  time.  For  these  meas¬ 
urements,  the  applied  voltage  was  9.4  kV,  the  streamer  cur¬ 
rent  probe  front  face  was  located  38  mm  from  the  capillary 
tip,  and  the  streamer  velocity  at  the  16  mm  location  was 
0.34  ±  0.05  mm/ns.  There  are  two  main  features  to  note  in 
Figure  7  that  may  help  elucidate  the  N2^  391  nm  and  N2 
337.1  nm  emission  intensity  behavior  for  the  helium  gas 
guided  streamer  discharges  with  air  or  N2-coflows.  First, 
note  that  both  the  rise  time  and  normalized  peak  emission 
intensities  of  the  N2^  391  nm  and  N2  337.1  nm  emission 
lines  are  time  delayed  by  5  ns  and  1 1  ns,  respectively,  rela¬ 
tive  to  the  peak  emission  intensity  of  the  He  587.6  nm  line. 
Based  on  the  streamer  velocity  measurement  for  the  applied 
voltage  condition,  and  assuming  that  the  peak  intensity  of 
the  He  587.6  nm  line  marks  the  streamer  head,  the  N2^ 
391  nm  and  N2  337.1  nm  peak  emission  intensities  follow 
approximately  2  mm  and  4  mm  behind  the  streamer  head, 
respectively.  Furthermore,  radial  slice  plots  of  the  He 
587.6  nm  emission  intensity  obtained  from  ICCD  imaging 
shown  in  Figures  8  and  9  indicate  that  the  streamer  head  at 
the  16  mm  location  is  always  located  well  within  the  He  gas 
core  flow  channel  since  the  optical  diameters  are  approxi¬ 
mately  one-third  of  the  capillary  tube  (2  mm)  inner  diameter. 
These  experimental  observations  indicate  that  the  head  of 
the  streamer  is  always  located  within  the  helium  gas  flow 
channel  and  co-flowing  nitrogen  plasma  emissions  are 
occurring  behind  the  streamer  head.  The  N2^  (B^Eu^) 


FIG.  7.  Spectrally  filtered  normalized  emission  intensity  traces  of  He 
587.6  nm,  He  706.5  nm,  N2^  391  nm,  and  N2  337.1  nm  transitions  obtained 
using  a  dual  fiber  coupled-PMT  optical  system  with  the  imaging  volume 
transverse  to  the  propagation  axis  and  located  16  mm  from  the  capillary  tip. 
The  inset  shows  normalized  ICCD  camera  radial  slice  plots  of  the  He 
587.6  nm  (solid  line)  and  N2^  391  nm  (dash-dot  line)  emission  intensities 
taken  under  similar  conditions  as  the  PMT  data. 


FIG.  8.  Normalized  radial  slice  plots  taken  from  ICCD  images  of  the  He 
587.6  nm  emission  intensity  at  gate  delay  times  when  the  streamer  head  was 
located  6  mm  (solid  line)  and  16  mm  (dashed  line)  from  the  capillary  tip. 
The  total  applied  voltage  was  6  kV  and  the  pressure  was  365  Torr. 

population  is  probably  produced  by  heavy  particle  energy 
transfer  from  helium  metastable  states  at  or  near  the  co¬ 
flowing  gas  interface  behind  the  streamer  head  because  the 
radial  distribution  of  the  reduced  electric  field  is  relatively 


FIG.  9.  Normalized  radial  slice  plots  taken  from  ICCD  camera  images  of  the 
He  587.6  nm  emission  intensity  at  a  gate  delay  time  when  the  streamer  head 
was  located  16  mm  from  the  capillary  tip.  In  both  figures,  the  line  shape 
FWHM  represents  the  streamer  optical  diameter.  In  Figure  9(a),  the  total 
applied  voltage  was  6kV  for  cell  pressures  170  Torr  (dashed  dotted  line), 
270  Torr  (dashed  line),  and  520  Torr  (solid  line).  In  Figure  9(b),  the  total 
applied  voltage  was  8kV  for  cell  pressures  270  Torr  (dashed  line),  520  Torr 
(solid  line),  and  760  Torr  (dash-dotted  line).  During  these  measurements, 
the  streamer  current  probe  was  located  23  mm  from  the  capillary  tip. 
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low  near  the  diffusional  mixing  layer  of  the  annular  co-flow 
N2  for  the  conditions  shown  in  Figure  7.  Thus,  it  is  unlikely 
that  both  the  core-flow  helium  and  annular  co-flow  nitrogen 
gas  will  be  simultaneously  ionized  at  the  streamer  head  for 
typical  helium  gas  guided  streamer  discharges.  We  note  that 
modeling  results  for  a  helium  gas  guided  streamer  having  a 
peak  axial  reduced  electric  held  of  70  Td  within  a  radially 
uniform  1%  air/helium  mixture  flow  channel  (see  Figure  16 
of  Ref.  33)  indicate  that  the  peak  electron  impact  excitation 
rate  for  N2  (C^Ilu)  v’=o  production  also  occurs  behind  the  He 
(^  =  3)  peak  excitation  rate.  The  second  feature  to  note  in 
Figure  7  are  the  N2  337.1  nm  and  N2^  391  nm  e-fold  decay 
times  relative  to  the  He  587.6  nm  peak  emission.  With 
regards  to  N2  337.1  nm  emission,  the  model  calculations 
shown  in  Figure  16  of  Ref.  33  indicate  that  both  the  He 
(«  =  3)  and  N2  (C^n  u)  State  excitation  rates  decrease  by  3 
orders  of  magnitude  within  0.05  mm  behind  their  peak  val¬ 
ues.  Assuming  that  both  of  these  electron  impact  excitation 
rates  turn  off  at  their  peaks  in  Figure  16  of  Ref.  33,  and  using 
the  effective  radiative  lifetime  of  3ns^^  for  the  N2(C^nu) 
state  quenched  solely  by  N2  (X^E^g),  we  estimate  that  the  e- 
fold  decay  of  the  N2  population  is  2.5  ns 

(^0.75  mm)  behind  the  He  (/2  =  3)  peak  excitation  rate. 
Under  similar  conditions  of  peak  electric  held,  our  data  in 
Figure  7  shows  that  the  e-fold  decay  time  relative  to  the  He 

587.6  nm  emission  peak  is  43  ns  (15  mm).  This  implies  that 
the  prolonged  decay  of  the  N2  337.1  emission  intensity  in 
Figure  7  is  due  to  excitation  through  heavy  particle  kinetics 
(dissociative  recombination  of  N4^).  We  note  that  the 
authors  of  Ref.  34  reach  the  same  conclusion  in  their  spec¬ 
trally  resolved  experimental  observations.  Similarly,  the  N2^ 
(B^Eu^)  effective  lifetime  is  approximately  0.2  ns  in  pure 
N2.^"^  In  Figure  7,  the  e-fold  decay  time  relative  to  the  He 

587.6  nm  emission  peak  is  19  ns  (6  mm),  implying  that  the 
N2^  391  nm  emission  within  the  channel  is  sustained  by 
energy  transfer  from  He  2  metastables. 

B.  Variable  pressure  data 

As  mentioned  in  Sec.  Ill  A,  the  ICCD  radial  and  axial 
emission  profiles  indicated  that  the  influence  of  the 
Laplacian  electric  held  extended  to  approximately  12  mm 
from  the  capillary  tip.  Figure  8  shows  normalized  radial  in¬ 
tensity  slice  plots  of  He  587.6  nm  emission  for  ICCD  gate 
delay  times  when  the  streamer  head  maximum  emission  in¬ 
tensity  was  located  6  mm  (solid  line)  and  16  mm  (dashed 
line)  from  the  capillary  tip  for  a  discharge  condition  of  6  kV 
applied  voltage  at  365  Torr.  Both  normalized  radial  intensity 
slice  plots  in  Figure  8  were  obtained  from  background  sub¬ 
tracted  ICCD  images  using  Image!  and  are  the  result  of  a  10 
pixel  binning  average  in  the  axial  direction  at  the  maximum 
emission  intensity.  At  6  mm,  the  intensity  profile  exhibits 
two  peak  emission  intensities  at  approximately  ±0.25  mm 
away  from  the  propagation  axial  centerline.  This  emission 
profile  indicates  that  the  streamer  has  a  ring-like  morphology 
that  is  formed,  initially,  along  the  inside  of  the  2  mm  diame¬ 
ter  dielectric  capillary  tube  wall  by  the  Laplacian  electric 
held  generated  at  the  ring  anode. As  the  streamer 
propagates  out  of  the  capillary  tip  to  6  mm  distance,  shown 
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in  Figure  8,  it  is  guided  by  the  ring  shaped  Laplacian  electric 
held  in  the  helium  jet  which  is  surrounded  by  the  N2  co-flow 
gas^^’^"^  and  driven  by  both  the  superimposed  Laplacian  and 
Poisson  electric  fields.  From  radial  and  axial  emission  inten¬ 
sity  slice  plots  obtained  from  ICCD  imaging,  we  observed 
that  the  two  off-axis  emission  intensity  peaks  converge  into 
a  single  peak  on  axis  occurring  approximately  10  mm  from 
the  capillary  tip.  However,  axial  intensity  slice  plots  of  the 
ICCD  emission  at  9kV  and  lOkV  applied  voltages  indicated 
that  the  Laplacian  electric  field  extended  up  to  approxi¬ 
mately  12  mm.  At  16  mm,  the  Laplacian  electric  field  is  neg¬ 
ligible  for  all  applied  voltage  and  pressure  conditions 
studied,  and  the  normalized  Gaussian-like  emission  intensity 
profile,  shown  in  Figure  8,  represents  emission  intensity  pro¬ 
duced  from  the  self-consistent  Poisson  electric  field  at  the 
streamer  head.  All  pressure  scaling  properties  of  the  streamer 
data  presented  in  this  section  have  been  studied  under  the 
conditions  of  Poisson  electric  field  dominated  propagation 
dynamics  located  16  mm  from  the  capillary  tip. 

Figure  9  shows  pressure  dependent  normalized  He 

587.6  nm  ICCD  emission  intensity  radial  slice  plots  obtained 
at  (a)  6kV  for  cell  pressures  170  Torr,  270  Torr,  and 
520  Torr,  and  (b)  8kV  for  cell  pressures  270  Torr,  520  Torr, 
and  760  Torr.  Each  radial  intensity  slice  plot  was  obtained  by 
binning  10  pixels  in  the  axial  direction  at  the  maximum 
emission  intensity  of  the  streamer.  The  optical  diameters  in 
Figures  9(a)  and  9(b)  have  been  denoted  by  the  FWHM. 
Note  that  the  optical  diameters  at  520  Torr  in  both  Figures 
9(a)  and  9(b)  correspond  to  those  plotted  in  Figures  5  and  6 
at  6  kV  and  8kV.  The  optical  diameter  measurements  as  a 
function  of  pressure  for  6kV  and  8kV  are  given  in  Table  11. 
At  6  kV,  the  optical  diameters  increase  slightly  (within  the 
13%  measurement  uncertainty)  and  can  be  considered  nearly 
constant.  At  8  kV,  the  increase  in  optical  diameters  is  approx¬ 
imately  14%  larger  compared  to  6  kV,  but  the  variation  is  still 
near  the  imaging  resolution  limit.  These  results  indicate  that 
the  optical  diameters  obtained  from  the  helium  guided 
streamer  are  essentially  invariant  within  the  pressure  range 
studied  at  6kV  and  8kV  applied  voltages.  Furthermore,  the 
streamer  velocities  are  also  invariant  to  optical  diameter,  simi¬ 
lar  to  that  of  Figure  6,  when  displayed  on  a  velocity-diameter 
plot  (not  shown).  The  invariance  of  the  optical  diameter  with 
pressure  for  the  helium  guided  streamer  again  exhibits  a 
strong  departure  from  the  predicted  inverse  gas  density 

TABLE  II.  Pressure  dependent  measurements  (see  Figure  9)  of  streamer  op¬ 
tical  diameters  and  estimates  of  peak  electron  number  densities  for  6  kV  and 
8kV  total  applied  voltages.  The  uncertainty  of  the  diameter  measurements 
are  ±13%  while  that  of  the  number  density  estimates  are  ±20%. 


Applied  voltage 

P  (Torr) 

6kV 

8kV 

D  (mm) 

ng  (x  10^%m“^) 

D  (mm) 

ng  (xl0^%m“^) 

170 

0.58 

8.3 

210 

0.60 

7.0 

270 

0.63 

6.7 

0.67 

8.2 

365 

0.63 

5.7 

0.70 

7.0 

520 

0.63 

5.5 

0.74 

6.0 

760 

0.75 

5.5 
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dependence  based  upon  similarity  scaling  and  confirmed  by 
measurements^"^’^^  and  modeling  in  both  N2  and  air. 
Furthermore,  the  invariance  of  optical  diameter  with  streamer 
velocity  is  also  in  strong  contrast  to  measurements and 
2-D  model  calculations^^  for  streamer  discharges  in  both  N2 
and  air. 

The  streamer  velocity  and  peak  streamer  current  pressure 
scaling  are  shown  in  Figure  10  as  a  function  of  inverse  pres¬ 
sure  for  both  6kV  and  8kV  applied  voltages  (squares  and 
circles,  respectively).  Note  that  the  data  shown  in  Figure  10 
correspond  to  the  same  data  set  presented  in  Figure  9  with  the 
addition  of  a  set  taken  at  6kV,  210Torr.  The  “best  ht”  poly¬ 
nomial  (see  Sec.  Ill  A)  that  describes  the  scaling  is  linear  with 
inverse  pressure  for  both  streamer  velocities  (solid  lines)  and 
peak  currents  (dashed  lines)  for  both  applied  voltages. 
Thus,  both  streamer  velocity  and  peak  streamer  current  vary 
inversely  proportional  to  pressure  at  6  kV  and  8  kV.  At  6  kV 
and  8kV  the  htted  slopes  for  the  velocity  are  [132.248 
±  17.129Torr*mm/ns]  and  [201.692  ±  19.656  Torr-mm/ns], 
respectively.  The  fitted  peak  current  slopes  at  6kV  and 
8kV  are  [5015.585  ±  717.576 Torr-mA]  and  [10240.298 
±  2071.687Torr*mA],  respectively.  This  inverse  pressure 
scaling  is  in  agreement  with  both  measurements  and  modeling 
performed  for  streamers  in  air.^^ 

IV.  ANALYSIS  AND  DISCUSSION 

In  this  work,  the  invariance  of  the  optical  diameter 
scaling  with  applied  voltage,  streamer  velocity,  and  gas 
pressure  exhibits  the  most  striking  difference  between  self- 
sustained  low  dielectric  strength  helium  guided  streamer 
with  nitrogen  co-flow  and  that  of  streamer  propagation  in 
homogeneous  air  or  nitrogen.  For  the  2  mm  I.D.  CDBD 
operating  with  2  SLM  helium  flow  at  520Torr,  peak 
streamer  current  measurements  with  applied  voltage  as 
shown  in  Figure  3  imply  that  diffusional  mixing  of  the  co¬ 
flowing  N2  gas  into  the  helium  core  flow  is  negligible  for 
self- sustained  streamer  propagation  up  to  approximately 
33  mm  from  the  capillary  tip.  For  each  repetitively  pulse 
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FIG.  10.  Streamer  velocities  (closed  symbols)  and  corresponding  peak  cur¬ 
rents  (open  symbols)  measured  as  a  function  of  inverse  pressure  for  total 
applied  voltages  of  6kV  (squares)  and  8  kV  (circles).  The  streamer  velocities 
were  measured  16  mm  from  the  capillary  tip  with  the  streamer  current  probe 
front  face  located  23  mm  from  the  capillary  tip.  The  solid  and  dashed  lines 
are  linear  best  fits  (see  text)  to  the  measured  speed  and  peak  current  data 
points,  respectively. 


discharge,  the  normalized  radial  slice  plots  shown  in 
Figure  9  indicate  that  the  Poisson  electric  field  driven  posi¬ 
tive  streamer  propagates  well  within  the  pure  helium  core 
flow  as  indicated  by  the  extent  of  the  optical  diameter  rela¬ 
tive  to  the  2  mm  capillary  tube  inner  diameter. 

Certainly,  the  higher  dielectric  breakdown  strength  of 
the  surrounding  N2  co-flow  confines  the  streamer  propaga¬ 
tion  to  within  the  helium  gas  channel  as  its  morphology 
transitions  from  ring-shaped  to  a  highly  localized  Gaussian- 
like  radial  profile. We  have  experimentally  verified  this 
by  substituting  N2  co-flowing  gas  with  helium  at  atmos¬ 
pheric  pressure,  with  the  result  that  the  helium  plasma 
extending  beyond  the  glass  capillary  is  no  longer  a  gas 
guided  streamer.  The  invariance  of  the  optical  diameter  ver¬ 
sus  applied  voltage,  streamer  speed,  and  gas  pressure  within 
the  pure  helium  core  flow  is  a  consequence  of  this  disparity 
in  the  gas  breakdown  strengths  but  more  importantly,  it  is 
due  to  relatively  higher  remnant  ionization  within  the  pure 
helium  core  compared  to  that  of  the  helium-nitrogen  mixing 
layer  (discussed  below).  If  the  timescale  of  electron-ion 
recombination  within  the  pure  helium  core  is  significantly 
longer  than  that  within  the  mixing  layer  during  the  pulse 
discharge  period,  then  successive  streamers  will  propagate 
through  remnant  pre-ionization  that  is  preferentially  higher 
within  the  pure  helium  core  relative  to  that  within  the  mix¬ 
ing  layer.  This  could  result  in  the  highly  localized  streamer 
propagation  within  the  pure  helium  core  as  shown  in  our 
experimental  results. 

We  have  constructed  a  0-D  onset  condition  model  by 
which  we  qualitatively  demonstrate  the  evolution  of  streamer 
localization  due  to  remnant  ionization  at  the  end  of  the  first 
discharge  pulse.  With  this  model,  electron  losses  are  com¬ 
pared  at  the  end  of  our  experimental  200  jus  pulse  repetition 
period  due  to  electron-ion  recombination  within  the  pure  he¬ 
lium  core  to  those  within  the  N2-He  mixing  region  (up  to  1  % 
N2  admixture)  for  the  experimental  condition  of  520Torr.  In 
Sec.  Ill  A,  the  estimated  electron  number  density  behind  the 
streamer  head  for  the  applied  voltage  range  6-lOkV  is 
approximately  6  x  lO^^cm”^  (see  Table  I)  corresponding  to 
a  maximum  electric  field  within  the  streamer  head  of  about 
llkV/cm  (65  Td).  Table  III  shows  the  elementary  reaction 
set  and  rate  coefficients  used  in  the  model  ’  “  for  the 
dominant  ionic  species  that  participate  in  (i)  direct  electron 
impact  ionization,  (ii)  ion  conversion,  and  (iii)  electron-ion 
recombination.  We  used  Bolsig+^^  to  calculate  the  direct 
ionization  rate  coefficients  (R1)-(R3)  for  65  Td,  520Torr 
conditions.  The  ionization  rate  coefficient  in  Table  I  for  he¬ 
lium  (R3)  is  an  average  of  the  ionization  coefficients  com¬ 
puted  for  pure  helium  and  a  1%  N2-He  admixture,  while 
those  for  N2^(X)  (Rl)  and  N2^(B)  (R2)  were  computed  for  a 
1%  N2-He  admixture. 

Before  discussing  the  numerical  calculations,  we  first 
considered  analytical  estimates  of  the  remnant  ionization 
(also  elucidating  the  dominant  recombination  mechanism)  at 
the  end  of  the  200  /is  pulse  repetition  period  for  the  limiting 
cases  of  (i)  pure  helium  and  (ii)  1%  N2-He  admixture  with 
background  gas  conditions  of  520Torr  and  300  K.  In  the 
case  of  pure  helium,  the  net  loss  rate  of  electrons  due  to 
recombination  is  (see  Table  III), 
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TABLE  III.  Elementary  reactions  and  rate  coefficients  used  in  the  zero-dimensional  model. 


Reaction  no. 

Reaction 

Rate  coefficient 

Reference 

Direct  electron  impact  ionization,  65  Td  and  520  Torr 

R1  e  +  N2^N2+(X)  +  2e 

1.76  X  10^’’  cm^/s  (1%  Nz/He) 

32 

R2 

e  +  N2^N2+(B)  +  2e 

1.04  X  10^  cm^s  (1%  Nz/He) 

32 

R3 

e  +  He  ^  He^  +  2e 

8.58  X  10“"  cm^/s  (avg.,  see  text) 

32 

Ion  conversion 

R4 

N2'^  +  2N2^N4++N2 

5  X  10“^^  cm®/s 

36 

R5 

He+  +  2He^He2'^  +  He 

6.5  X  10“^^  cm'’/s 

37 

R6 

He2+  +  N2^N2+  +  He2 

1.4  X  10^^  cm^/s 

39 

Electron-ion  recombination 

R7 

e  +  N2+^N  +  N 

2.8  X  10^^(300  K/Te)°-^  cm^/s 

36 

R8 

e  +  Nj'^^N  +  NC^D) 

2  X  10“’  (SOOK/Te)'’-^  cm’/s 

36 

R9 

e  +  Hej'^^He  +  He 

3.3  X  10““’cmVsat300K 

38 

RlO 

e  +  N/^Nz  +  Nz 

2  X  10“'^  (SOOK/Te)”'^  cm^/s 

37  and  43 

Rll 

2e  +  He2+^He'"  +  He  +  e 

4  X  10“’"  (Tg/Te)”*  cm"/s 

37 

R12 

e  +  He  +  Hez^  ^  He”  +  2He 

5  X  10“”  (Tg/Te)  cm®/s 

37 

R13 

2e  “t“  He"*"  — ^  He  “t“  e 

8  X  10“’"  (Tg/ry  cm"/s 

37 

R14 

2e  +  N2'^^e  +  N2 

1  X  10“'"  (300K/Te)"''^  cm"/s 

36 

R15 

e  +  Nz'^  +  Nz^Nz  +  Nz 

6  X  10“”  (300K/Te)''^  cm"/s 

36 

40 


d/? 

=  — ^9^//e+  —  l^n^He+  ~  l^l2^He^He+^e  “ 

(2) 

where  13^  are  the  recombination  rate  coefficients  (R9), 
(R11)-(R13)  in  Table  III  and  represent  the  ion,  electron, 
and  neutral  gas  number  densities.  Note  that  the  electron  den¬ 
sity  (6  X  lO^^cm”^)  remains  constant  during  ion  conversion 
so  it  then  becomes  the  initial  condition  for  the  recombination 
process.  In  addition,  we  have  assumed  that  the  electron  tem¬ 
perature  during  recombination  is  equal  to  the  neutral  gas 
temperature,  300  K.  The  timescale  for  ion  conversion  is  typi 
cally  much  faster  than  that  of  electron-ion  recombination. 
At  520Torr,  He^^He2^  (R5)  ion  conversion  timescale  is 
approximately  56  ns  while  that  of  the  atomic  ion  recombina¬ 
tion  term  (R13)  is  approximately  50  /is  for  an  initial  electron 
number  density  n^  o  =  6  x  10^^  cm“^.  Thus,  all  He^  ions  con¬ 
vert  to  He2^  prior  to  recombination  so  that  the  last  rate  term 
in  Eq.  (2)  is  negligible.  Charge  conservation  requires  that  the 
He2^  density  equals  that  of  the  electrons  during  recombina¬ 
tion,  allowing  us  to  estimate  the  magnitudes  of  the  three 
remaining  terms  in  Eq.  (2).  The  rate  of  three-body  dissocia¬ 
tive  recombination  (R12)  dominates  those  of  (R9)  and  (Rll) 
by  at  least  an  order  of  magnitude,  so  that  Eq.  (2)  has  the  ap 
proximate  solution 


40 


^e(0 


^e,0 


1  +^e,0i5l2^He^’ 


(3) 


where  =  0)  =  ^e,o  =  6  x  10^^  cm“^  and  riue  is  the  neutral 
background  helium  density.  Note  that  the  initial  time  can  be 
defined  as  ^  =  0  ns  with  good  approximation  because  of  the 
comparatively  long  electron-ion  recombination  timescale. 
Evaluating  Eq.  (3)  just  before  the  next  excitation  pulse  at 


200  JUS,  we  find  that  ^  5  x  10  cm“  within  the  pure  he¬ 
lium  core. 

The  net  electron  recombination  rate  equation  for  a  1% 
N2-He  mixture  at  520Torr,  300  K  (reactions  R7-R15)  is 


d^e 

“dT 


-  -  ^12«He«He+«e  “  ^13«He+«e 

—  ^14«N+^e  “  ^15^N2%tWe, 


(4) 


where  (3^^^  =  l3'j-\-  (3^.  Equation  (4)  can  be  greatly  simplified 
because,  again,  the  ion  conversion  and  electron-ion  recombi¬ 
nation  occur  nearly  in  series  on  very  different  timescales. 
The  ion  conversion  pathway  in  this  case  is  He^^He2^ 
^  N2^  ^  N4^.  Eor  these  conditions,  the  conversion  of 


He2^  ^  N2^  is  much  faster  (4  ns)  than  that  of  N2^  ^  N4^ 
(1  /is).  By  the  end  of  the  ion  conversion  process,  N4^ 
becomes  the  dominant  ion  and  is  equal  to  the  electron  den¬ 
sity  during  recombination.  We  note  that  total  ion  conversion 
to  N4^  occurs  on  a  similar  timescale  as  that  of  the  N4^  disso¬ 
ciative  recombination  reaction  (RIO),  but  the  numerical  solu¬ 
tion  discussed  below  indicates  that  this  is  a  good 
approximation.  The  net  electron  number  density  rate  equa¬ 
tion  has  the  solution 


ne(0  = 


^e,0 


1  +  ^e,0^lo^ 


(5) 


where  n^it  =  0)  =  n^  o  =  6  x  lO^^cm”^.  Evaluating  Eq.  (5)  at 
the  end  of  the  pulse  period,  we  find  ^  2.5  x  lO^cm”^ 
within  the  1%  N2-He  mixing  layer.  These  analytic  estimates 
for  electron-ion  recombination  at  520Torr  indicate  that  the 
remnant  ionization  within  the  pure  helium  core  at  the  end  of 
the  200  pis  pulse  period  is  an  order  of  magnitude  greater 
than  that  within  the  1%  N2-He  mixing  layer. 

Eigures  11  (a)- 11(c)  show  the  results  of  solving  the  set  of 
time  dependent  net  rate  equations  for  electron,  He^,  He2^, 
N2^,  and  N4^  densities  at  520Torr,  300  K  neutral  gas  temper¬ 
ature,  and  pre-ionization  density  of  lO^cm”^  for  (a)  pure  he¬ 
lium,  (b)  0.01%  N2-He,  and  (c)  1%  N2-He  admixtures  using 
the  4th  Order  Runge-Kutta  method."^^  In  this  0-D  model,  we 
have  assumed  that  the  electron  and  gas  temperatures  are  equal 
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FIG.  11.  Zero-dimensional  modeling  results  of  electron  (solid  line),  He^ 
(dash),  He2^  (dash-dot),  N2^  (dotted),  and  N4^  (dash-dot-dot)  number  den¬ 
sities  vs.  time  for  the  cases  of  (a)  pure  helium,  (b)  0.01%  N2-He  admixture, 
and  (c)  1%  N2-He  admixture. 

during  electron-ion  recombination.  We  have  chosen  the  value 
of  0.01%  N2-He  in  Figure  11(b)  because  it  approximately 
demarcates  the  admixture  for  which  the  nitrogen  ion  species 
begin  to  play  a  significant  role  in  ion  conversion  and  recombi¬ 
nation  and  therefore,  becomes  an  upper  bound  for  which  the 
core  flow  can  be  considered  as  “pure”  helium.  In  each  case 
shown  in  Figures  ll(a)-ll(c),  we  have  modeled  the  passing 
of  the  self-propagating  electric  field  at  65  Td  as  a  square  pulse 
during  which  we  have  solved  the  rate  equations  (R1)-(R6) 
until  the  electron  number  density  reached  6  x  lO^^cm”^.  The 
time  to  reach  this  electron  density  (solid  line)  was  approxi¬ 
mately  5  ns  and  was  practically  independent  of  nitrogen 
admixture  (note  that  this  computed  rise  time  varied  between 
approximately  2  and  10  ns  for  pre-ionization  densities 
between  10^^  and  10"^cm“^,  respectively).  This  time  is  much 
shorter  than  the  measured  streamer  current  rise  time  of  70  ns 
shown  in  Figure  2  (Sec.  Ill  A).  Since  the  streamer  current  rise 
time  may  be  impacted  due  to  electron  drift  through  the 
streamer  channel  and  Pearson  probe  circuit,  we  have  used  the 
results  of  the  He  (3^D— 2^P)  587.6  nm  emission  intensity  rise 
time  acquired  using  the  side-view  PMT  shown  in  Figure  7  as 


well  as  ICCD  camera  axial  slice  plots  (not  shown)  taken  with 
a  5  ns  gate  for  comparison.  The  He  587.6  nm  rise  time  using 
the  PMT  in  Figure  7  was  about  3  ns  while  that  using  the 
ICCD  camera  was  5  ns.  Hence,  Figures  ll(a)-ll(c)  show  that 
the  calculated  time  to  reach  an  electron  number  density  of 
6xl0^^cm“^,  which  was  estimated  using  experimentally 
measured  quantities  expressed  by  Eq.  (1),  are  in  qualitative 
agreement  with  that  of  two  independent  He  (3^D— 2^P) 
587.6  nm  emission  intensity  rise  time  measurements. 

Figure  11(a)  shows  the  results  for  the  case  of  pure  he¬ 
lium  in  which  all  of  the  He^  population  (dashed  line) 
converted  to  He2^  (dash-dotted  line)  after  1  jus,  and  the  rem¬ 
nant  ionization  is  5  x  lO^^cm”^  after  2  x  10“^  s  (200  jus),  in 
agreement  with  the  approximation  using  Eq.  (3).  Eor  pure 
helium,  only  the  three-body  electron-ion  recombination  reac¬ 
tion  (R12)  had  any  effect  on  the  electron  decay  rate,  which  is 
again  consistent  with  the  approximation  for  obtaining  Eq. 
(3).  Eor  the  nitrogen  admixtures  shown  in  Eigures  11(b)  and 
1 1(c),  we  did  not  use  reaction  (R12)  because  helium  metasta¬ 
bles  are  rapidly  quenched  through  Penning  ionization  with 
subsequent  nitrogen  ion  conversion  and  dissociative  recom¬ 
bination.  Eigure  11(b)  shows  the  results  for  0.01%  N2-He 
admixture.  In  this  case,  both  He^  and  He2^  are  converted 
into  N2^  (dotted  line),  which  has  a  maximum  of 
4  X  10^^  cm“^  at  about  1  jus.  Note  that  the  N4^  (dash-dot-dot¬ 
ted  line)  population,  produced  via  ion  conversion  reaction 
(R4),  is  less  than  lO^cm”^  and  that  the  dissociative  recombi¬ 
nation  reaction  R(10)  has  no  effect  on  the  electron  decay 
rate.  The  electron  decay  rate  is  due  solely  to  N2^  dissociative 
recombination  reactions  (R7)  and  (R8)  and  the  remnant  ioni¬ 
zation  at  200  /is  is  1  X  10^^cm“^  (five  times  smaller  than  that 
of  pure  helium).  Eigure  11(c)  shows  the  results  for  the  1% 
N2-He  admixture.  In  this  case,  the  concentration  of  neutral 
nitrogen  is  sufficient  that  the  sequential  ion  conversion  path¬ 
way  He2^  ^  N2^  leads  to  N4^  (R4)  becoming  the  dominant 
ion  (2  X  lO^^cm”^)  within  1  ^s.  The  dissociative  recombina¬ 
tion  of  N2^  (R7  and  R8)  was  found  to  affect  the  onset  of 
electron  decay  at  about  0. 1  /is,  but,  as  the  dominant  ion  spe¬ 
cies  concentration  switched  from  N2^  to  N4^  at  approxi¬ 
mately  1  /IS,  dissociative  recombination  (RIO)  became  the 
dominant  electron  loss  mechanism.  The  remnant  ionization 
in  Eigure  11(c)  at  200 /is  is  2.5  x  lO^cm”^  and  is  the  same 
result  obtained  using  the  approximate  Eq.  (5). 

V.  CONCLUSION 

The  0-D  electron-ion  recombination  model  shows  that 
within  the  gas  mixing  layer  where  the  nitrogen  concentration 
is  1%,  the  remnant  ionization  after  200 /is  is  more  than  an 
order  of  magnitude  smaller  than  that  within  the  pure  helium 
channel,  assuming  that  the  reduced  electric  field  at  the 
streamer  head  is  65  Td  and  the  electron  number  density 
behind  the  streamer  is  on  the  order  of  6  x  lO^^cm”^.  If  the 
self-propagating  electric  field  at  the  streamer  head  has  a  sim¬ 
ilar  Gaussian-like  profile  as  shown  by  the  He  (3^D— 2^P) 
587.6  nm  emission  measurements  in  Eigure  9,  then  the  rem¬ 
nant  ionization  in  the  mixing  layer  will  be  lower  than  these 
estimates  compared  to  the  pure  helium  channel.  This  sug¬ 
gests  that  under  our  experimental  conditions  of  5  kHz  pulse 
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repetition  rate,  16  mm  from  the  capillary  tip,  it  is  the 
electron-ion  recombination  rates,  rather  than  the  direct  ioni¬ 
zation  rates  within  the  1%  mixing  layer, that  maintain  a 
dominant  conductive  channel  within  the  pure  helium  gas 
core  which  guides  streamer  propagation  and  prohibits 
streamer  branching  with  increase  in  applied  voltage.  In  other 
words,  for  typical  CDBD  reduced  electric  helds,  the  lower 
dielectric  breakdown  strength  of  the  helium  gas  channel 
compared  to  the  surrounding  nitrogen  annular  co-flow  and 
the  electron  loss  rates  (rather  than  direct  ionization  produc¬ 
tion  rates)  become  the  proper  boundary  conditions  that  lead 
to  the  observed  invariance  of  the  optical  diameter  with 
applied  voltage,  streamer  velocity,  and  gas  pressure. 

Without  diffusional  mixing  of  air  or  nitrogen,  guided  he¬ 
lium  streamer  model  simulations^^  indicate  that  the  initial 
externally  applied  Laplacian  electric  field  is  responsible  for 
the  formation  and  guiding  of  the  ring-shaped  streamer  into 
the  helium  core  flow  outside  the  dielectric  tube,  rather  than 
direct  electron  impact  ionization  at  the  1  %  helium- air  mixing 
layer. The  simulations^^  also  show  that  the  high  off-axis 
ionization  region  converges  towards  the  axis  faster  in  time  as 
the  initial  uniform  pre-ionization  density  decreases  from 
lO^cm”^  to  10"^cm“^,  which  is  consistent  with  a  larger  on- 
axis  electric  field.  With  the  effect  of  diffusional  mixing, 
the  confinement  of  the  ring-shaped  ionization  region  may 
also  be  a  consequence  of  electron-ion  recombination,  rather 
than  direct  electron  impact  ionization,  in  addition  to  the  dif¬ 
ference  in  breakdown  strengths  of  helium  versus  nitrogen. 
The  Laplacian  electric  field  is  then  responsible  for  guiding 
the  streamer  into  the  helium  core  flow  upon  exiting  the 
dielectric  tube  rather  than  the  1%  helium-nitrogen  mixing 
layer. 

Finally,  the  temporal  variation  of  both  rise  times  and 
decay  times  of  emission  intensity  plots  in  Figure  7  show 
that  only  helium  emission  lines  are  excited  by  direct  elec¬ 
tron  impact.  Both  N2  C-B  and  N2^  B-X  emission  lines  are 
produced  by  heavy  particle  kinetics.  The  N2^  B-X  emission 
line  is  most  likely  produced  by  the  energy  transfer  from 
He(2  ^S)  +  N2  ^  N2^  B  +  He  with  measured  rate  constant 
1  X  10“^  cm^/s."^^  The  N2  C  state  excitation  is  produced  by 
N4+  +  e^  N2  (C^n  u)  +  N2  with  measured  dissociative 
recombination  rate  constant  of  2.6  x  10“^  cm^/s  at  300 
The  N2  C  ^  B  337.1  nm  emission  marks  the  dissociative 
recombination  reaction  of  N4^  ions  in  streamer  channel 
behind  the  streamer  head. 
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